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Abstract 



C^ , We analyze the effect of the simultaneous presence of correlated additive 

a, 

I , and multiplicative noises on the stochastic resonance response of a modulated 

o 
o 



bistable system. We find that when the correlation parameter is also mod- 
ulated, the system's response, measured through the output signal-to-noise 



^^ I ratio, becomes largely independent of the additive noise intensity. 
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The phenomenon of stochastic resonance (SR) has attracted considerable interest in the 
last decade due, among other aspects, to its potential technological applications for optimiz- 
ing the output signal-to-noise ratio (SNR) in nonlinear dynamical systems. The phenomenon 
shows the counterintuitive role played by noise in nonlinear systems as it contributes to en- 
hance the response of a system subject to a weak external signal. There is a wealth of 
papers, conference proceedings and reviews on this subject, Ref. |jl| being the most recent 
one, showing the large number of applications in science and technology, ranging from pa- 
leo climatology, to electronic circuits, lasers, and noise-induced information flow in sensory 
neurons in living systems. 

A tendency shown in recent papers, and determined by the possible technological appli- 
cations, points towards achieving an enhancement of the system response (that is obtaining 
a larger output SNR) by means of the coupling of several SR units 0-0 in what conforms an 
"extended medium" |^. Another aspect that has also attracted interest is finding a system 
where the SNR becomes largely independent of external parameters such as noise intensity 

ii- 

In this work we present some results related to the latter aspect. In order to reach 

such a goal we study a bistable system subject to both an additive and a multiplicative 

noise but, at variance with the work in Ref. [§], we consider that both noises are correlated. 

In addition to the modulation of the bistable potential by a weak external signal, we also 

consider the effect of a modulation of the correlation between both noises. We show that 

the result of this contribution is to widen the maximum of the SNR as a function of the 

additive noise intensity, making the detection of the signal less sensitive to the actual value 

of that noise. It is worth remarking here that on one hand, the additive (external) noise 

can be (naturally) assumed as white, while on the other hand, the multiplicative (internal) 

one, that generally involves characteristic time scales of the system, would not necessarilly 



be white. However, as indicated in |T^, as a first step in this study, we can approximate the 
multiplicative colored noise by a white one. Here we exploit the results of Madureira et al. 
||T0| where it was shown that the activation rate in a bistable system subject to correlated 



additive and multiplicative noises is dramatically suppressed or enhanced according to the 
sign of the correlation. The model equation we consider corresponds to a onedimensional 
bistable system described by the Langevin equation 

x = e{t)+x-x^ + x^i{t)+^2{t), (1) 

where e(t) = Eq cos{Qt), Eq and Q are the intensity and frequency of the potential modulation 
respectively. The additive and multiplicative Gaussian white noises, indicated by .^i(t) and 
^2(t), satisfy 

(6(t)) = 0, z = l,2 

{Ut)Ut')) = 26{t-t')Q 
{UtMt')) = 26{t-t')D 



{UtMt')) = 2p{t)6{t-t')^DQ. (2) 

The correlation intensity between both noises is indicated by the parameter p{t), that fulfills 
the condition \p\ < 1. The associated the Fokker-Planck equation (Stratonovich prescrip- 
tion) 



dtP{x,t)= -d^ ix-x'^ + Qx + p{t)^QD-E{t)] P{x,t) 

+ Ddl^[(l + Rx'^ + 2p{t)VR)P{x,t)]. (3) 

In what follows we assume that p is a time dependent function having the form p{t) = 
Pq cos{Qpt). Here, and in order to simplify the analysis, we consider the case when Qp = Q, 
the most general case will be discussed elsewhere |jlT[ . 



In order to perfom the evaluation of the correlation function and power spectral density 
needed to obtain the SNR, we exploit the results of the two-state model p!^JT3[] . Those 
authors have reduced the problem of obtaining the SNR of a nonlinear and essentially 
bistable system subject to a weak periodic signal, to a description where the transitions 
occur between the two minima of the deterministic potential. The main approximation is 
an adiabatic-like one, corresponding to the assumption that the relaxation times around each 



minima are much shorter than any other characteristic time, such as the Kramers time for 
transitions between the two stable points or the inverse of the signal frequency Q. It is also 
important to note here that, if the modulation is done around p 7^ 0, it becomes necessary 
to extent the two state approach in order to include the asymmetry of the potential. 

In the absence of any signal, the deterministic potential of the system has two minima 
located at the points x± = ±1. It has been shown that such minima are not coincident with 



the maxima of the stationary probability distribution [|1J . 

Let us call W+{t) and W-{t) the nonstationary transition rates from the state x+ to 
X- and from the state x_ to x+ respectively. As indicated before, the signal is introduced 
through a simultaneous modulation in the potential and in the correlation. We may note 
here that the last contribution does not modify the deterministic potential barrier. However, 



what we find is that the rates W± do change, i.e. in |jTO| it was shown that if p > the rate 
Wj^ increases while W- decreases, and viceversa. 

If both modulations are small in comparison with the barrier height, that is po ^ 1 and 
£0 ^ V{0) — V{±1), then it is possible to make a Taylor expansion of the functions W±{t). 
We thus obtain 

Wi{t) = l{WoT («p + «e) cos(fit) + 0{sl) + 0{pI) ■■■), (4) 

where Wo is the transition rate evaluated in absence of modulation, po = Eq = 0. The terms 
ttp and as are given by 

as dWo 



2 de 

as dWo 



So (5) 



e=0 



Po- (6) 

p=0 



2 dp 

The transition rate Wo could be obtained calculating the mean-first-passage-time T{R,p) 
n| , yielding 



^ = T{R,p) = ^£ dxH{x) exp Mx)/D] 



r dyH{y)exp[-^y)/D]. (7) 

J —cx> 



where R = Q/D, and the function H{x) is 

H(x) -- 



1 



V 1 + i?a;2 + 2pVRx 

The determination of the effective potential $(x) is made studying the stationary probabihty 
density, which is given by 



$(a;) = r H{x'f (l + Rx''^ + (^o cos(fit) + 2p/R) x') dx' . 

We obtained the derivatives dWo{R,p)/de and dWo{R, p)/dp from the evaluation of 
dT{R,p)/de and dT{R,p)/dp respectively, through 

dWoiR,p) 1 dT{R,p) 



de T{R,py de 

dWo{R,p) 1 dT{R,p) 



(8) 
(9) 



(ip T{R,py dp 

The results for the derivatives of T{R,p), after a long but straightforward calculation, are 
,dT 



D' 



de 



dxHix)eic£>\^{x)lD] T dyHiy) exp [-$(|/)/L)] T dyHiyfy 
-1 J-i J-i 



+ f^ rfy/7(t/)exp[-<l>(y)/D]y'' dzH{zfz, 



(10) 



for the contribution of the potential modulation, while for the contribution of the correlation 
modulation we obtain 

. — = - / (ia;iJ(a;)exp [$(a;)/D] X iJ(a;)^a; / dyH(y\&^^\-^(y)lD\ 
vQ-D dp J-i J-i 

+ ^ £ dyHiy) exp [-Hy)/D] £ dyy'il - y^)H{yf 

-£ dyH{yfy + ^f^ dy H{y)exp[-^y)/D] f^z'{l - z')H{z)\ (11) 

Once «£ and ap are determined, we can obtain the auto-correlation function, which is 
given by 



{x{t)\x{t + r)) 



-Wo\t\ 



1- 



(«£ + dp)^ cos(f2t — 
a + n 



+ 



[a, + apf (cos(fir) + cos [fi(2t + r) + 20])) 
5 



(12) 



From {x{t)\x{t + r)) the power-spectrum density (psd) {S^u)) is readily obtained as 

2Wo 



(5H)* 



^ _ {ap + aJ2 



W^ + o;2 



7T[ap + a^y 



+ 2(i^^^--^) (1^) 



where we have written only the expression for positive values of uj ||T^ . 
Then, to determine the output SNR TZ, we use the standard definition 

where 5"^ and Sn are the (integrated) psd with and without signal respectively, evaluated 
at the modulation frequency. Here, the parameter A is introduced in order to tune the 
theoretical result when compared with a numerical simulation or an experiment, A being 
related to the bandwidth of the sampling frequencies. It is clear that the evaluation of TZ is 
critically dependent on the evaluation of dWo/de and dWo/dp. For this reason the results 
of Eqs. (p!0| , p!T]) have been numerically tested for a large range of values of the different 
parameters. 

One of the goals of the present work was to look for a larger independence of the SNR 
from the external parameters, particularly the additive noise intensity. In Fig. 1 we show 
the results for 7?. as a function of the noise intensities D and Q, for a fixed value of po and 
different values of eo- The widening of the peak as a function of D for a region of values 
of Q is apparent. However, this widening occurs at the expense of a (small) reduction in 



the maximum of TZ. The limiting case for Q ^ 0, reducing to the usual SR behaviour [|13 
should be compared with the case when Q 7^ to see clearly the widening effect. Also, the 
behaviour for D — > (and Q y^ 0) reduces to the case studied in 0. Such a widening of TZ 
is due to the additive dependence of the ouput signal on Pq. 

Another novel aspect of our results is that when Eq —>■ 0, we have a SR phenomenon 
purely due to the modulation of the correlation parameter that is more localized in the 
{D, Q) parameter space. When one or both of these parameters goes to or 00 we find that 

The results obtained by modifying only the intensity of correlation modulation amplitude 



are depicted in Fig. 2. In this case, the results are essentially the same as in the previous 
figure. 

It is worth to note here that in Figs. lb,c and 2b, keeping constant D, the SNR grows 
as Q increases. This remarkable fact implies that the additive SR (the case (5 = 0) could 
be enhaced by adding a multiplicative noise. 

We have also studied the eff'ect of changing the frequency modulation Vt on SNR, observ- 
ing the same dependence on this parameter, as in other SR sytems. For very low frequen- 
cies the effect is slightly larger than for higher frequencies, but it is essentially frequency- 
independent. 

Finally, in Fig. 3, we depict a contour plot of the SNR as a function of po and ^o- 
We fixed D^ Q and the frequency modulation Vt. A paraboloid-like dependence on those 
variables is obtained, with TZ increasing as any of the variables increases too. 

In conclusion, we have shown that the simultaneous modulation of the potential and 
the correlation parameter between the additive and multiplicative noises induces a widening 
of the output SNR, making the system response less dependent on the precise value of the 
additive noise. It is worth noting that if the modulation frequencies for the potential and the 
correlation are too different, what we find is a superposition of the SR effect associated to 
each separated modulation, that is: two differenciated peaks in the SNR at the corresponding 
frequencies. 

The present result could be of relevance for both technological and biological systems, 
as in electronic signal detectors or sensory systems one wishes that the detection capacity 
shall be as little dependent as possible on the (usually external) additive noise, the (usually 
internal) multiplicative noise being a kind of tuning parameter. A more complete study of 
the dependence of this effect on the modulation frequency of the correlation will be presented 
elsewhere |Tl[ . 

Acknowledgments: The authors thank P. Hanggi and G. Nicolis for useful comments and 
V. Grunfeld for a revision of the manuscript. Partial support from CONICET, Argentina, 
through grant PIP Nro.4593/96, ANPCyT, Argentina, through grant 03-00000-00988, and 



CEB, Bariloche, Argentina, are also acknowledged. 



REFERENCES 

[1] L. Ganimaitoni, P. Hanggi, P. Jung and F. Marchesoni, Rev. Mod. Phys. 70, 223 (1998). 

[2] A. Bulsara and G. Schmera, Phys. Rev. E 47, 3734 (1993); P. Jung, U. Behn, E. 
Pantazelou and F. Moss, Phys. Rev. A 46, R1709 (1992). 

[3] P. Jung and G. Mayer-Kress, Phys. Rev. Lett. 74, 208 (1995). 

[4] J. Lindner, B. Meadows, W. Ditto, M. Inchiosa and A. Bulsara, Phys. Rev. Lett. 75, 3 
(1995); ibid, Phys. Rev. E 53, 2081 (1996). 

[5] H. S. Wio, Phys. Rev. E 54, R3075 (1996). 

[6] F. Marchesoni, L. Ganimaitoni and A. Bulsara, Phys. Rev. Lett. 76, 2609 (1996). 

[7] H. S. Wio and F. Castelpoggi, Proc. Conf. UPoN'96, C.R.Doering, L.B.Kiss and 
M.Schlesinger Eds. (World Scientific, Singapore, 1997), pg. 229; F. Castelpoggi and 
H.S. Wio, Europhysics Letters 38, 91 (1997); F. Castelpoggi and H. S. Wio, Phys. Rev. 
E 57, 5112(1998); M. Kuperman, H.S. Wio, G. Iziis and R. Deza, Phys. Rev. E 57, 
5122 (1998); M. Kuperman, H.S. Wio, G. Iziis, R. Deza and F. Castelpoggi, Physica A 
(1998) in press. 

[8] J.J. Colhns, C.C. Chow and T.T. Imhoff, Phys. Rev. E 52, R3321 (1995); J.J. Collins, 
C.C. Chow, A.C. Capela and T.T. Imhoff, Phys. Rev. E 54, 5575 (1996). 

[9] L. Gammaitoni, F. Marchesoni, E. Menichella-Saetta and S. Santucci, Phys. Rev. E 49, 
4878 (1994). 

[10] A.J.R. Madureira, P. Hanggi, H.S. Wio. Phys. Lett. A 217, 248 (1996). 

[11] C.J. Tessone and H.S. Wio, in preparation. 

[12] C. Nicolis, Tellus 34, 1 (1982). 

[13] B. McNamara, K. Wiesenfeld, Phys. Rev. A 39, 4854 (1989). 



[14] A. Fulinski and T. Telejko, Phys. Lett. A 152, 11 (1991). 



10 



FIGURES 
FIG. 1. Here we depict the SNR as function of the additive and multiplicative noise intensities 

D and Q. On the left side we show the 3D plots associated with the contour plots on the right 

side. We have fixed po = 0.05 and the frequency modulation 0, = 0.005, while the different cases 

correspond to: (a) Eq = 0.01, (b) = 0.05, (c) = 0.10. For the sake of clarity, the point of view in 

the left side 3D plots, is from the upper left corner in the contour plots. 

FIG. 2. We show the SNR in {D,Q) space. Once again, the left side 3D plots correspond to 
the contour plots on the right side, and the point of view for the 3D plots corresponds to the upper 
left corner in the contour plots. We fixed the frequency modulation 0, = 0.005, the amplitude of 
potential modulation eq = 0.05, and varied the potential amplitude: (a) corresponds to po = 0.01; 
(b) = 0.10. 

FIG. 3. Contour of 7^ as a function of the modulation intensities po and Eq. The values for the 
noises are D = 0.25, Q = 0.10 and for the frequency modulation Q = 0.005. 
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